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Abstract 

A systematic study is performed of fully integrated particle multiplicities in 
central Au-Au and Pb-Pb collisions at beam momenta of 1.7 A GeV/c, 11.6^4 GeV/c 
(Au-Au) and 158A GeV/c (Pb-Pb) by using a statistical-thermal model. The close 
similarity of the colliding systems makes it possible to study heavy ion collisions 
under definite initial conditions over a range of centre-of-mass energies covering 
more than one order of magnitude. In order to further study the behaviour of 
strangeness production, an updated study of Si-Au collisions at 14. 6A GeV is also 
presented. The data analysis has been performed with two completely independent 
numerical algorithms giving closely consistent results. We conclude that a thermal 
model description of particle multiplicities, with additional strangeness suppression, 
is possible for each energy. The degree of chemical equilibrium of strange particles 
and the relative production of strange quarks with respect to u and d quarks are 
higher than in e + e~, pp and pp collisions at comparable and even at lower energies. 
The behaviour of strangeness production as a function of centre-of-mass energy 
and colliding system is presented and discussed. The average energy per hadron in 
the comoving frame is close to 1 GeV per hadron despite the fact that the energy 
increases more than 10-fold. 

PACS:24.10.Pa,25.75.Dw,25.75.-q 



1 Introduction 



After scouring results from relativistic heavy ion collisions at many different energies over 
several years some common traits are starting to emerge. Indeed, statistical-thermal models 
have proved to be able to reproduce relative particle multiplicities in a satisfactory manner 
by using two or three relevant parameters: temperature, baryon chemical potential and a 
possible strange-quark suppression parameter, 7 S 0. Such an analysis has been performed 
by many authors for heavy ion collisions data from CERN SPS, from Brookhaven AGS and 
also from GSI SIS. In this paper we present a simultaneous analysis of data from several 
different collisions, with emphasis on the similarity of the colliding system in order to study the 
behaviour of parameters as a function of centre-of-mass energy within one framework. Hence, 
we have focussed our attention on central Au-Au collisions at beam momenta of 1.74 GeV/c 
(SIS) §, 11.64 GeV/c (AGS) [§ and on central Pb-Pb collisions at 1584 GeV/c (SPS) §. 
As far as the choice of data (and, consequently, colliding system) is concerned, our leading 
rule is the availability of full phase space integrated multiplicity measurements because a pure 
statistical-thermal model analysis of particle yields, without any consideration of dynamical 
effects, may apply only in this case ||. Such data, however, exist only in a few cases and 
whenever legitimate we have extrapolated spectra measured in a limited rapidity window to full 
phase space. The use of extrapolations is more correct than using data over limited intervals of 
rapidity, especially in the framework of a purely statistical-thermal analysis without a dynamical 
model. Moreover, the usual requirement of zero strangeness (S = 0) demands fully integrated 
multiplicities because strangeness does not need to vanish in a limited region of phase space. 

A point of considerable interest in heavy ion collisions is the enhanced production of strange 
quarks per u, d quark with respect to elementary collisions [[]] like e + e~, pp , pp . This could 
be related to properties of the system at the parton level prior to hadronisation |l|, ||, [7|, |8j . In 
order to further study strangeness production and enhancement at low energy, we also present 
a new analysis of Si-Au collisions at 14.64 GeV (AGS) using only multiplicities obtained from 
fully integrated phase space distributions. This also allows to cross-check results of previous 



analyses [pi nil 13 1 performed using limited rapidity interval data. In particular, we have 



included the 4ir pion multiplicity |TTJ and results presented in [JL2]. In order to assess the 



consistency of the results obtained, we have performed the statistical-thermal model analysis 
by using two completely independent numerical algorithms whose outcomes turned out to be 
in close agreement throughout. 



Similar analyses have been recently made by other authors (see e.g. [13, |14j); however, both 
the model and the used data set differ in several important details, such as the assumption 
of full or partial equilibrium for some quark flavours, the number of included resonances, the 
treatment of resonance widths, inclusion or not of excluded volume corrections, treatment of 
flow, corrections due to limited rapidity windows etc. Because of these differences it is difficult to 
trace the origin of discrepancies between different results. We hope that the present analysis, 
covering a wide range of beam energies using a consistent treatment, will make it easier to 
appreciate the energy dependence of the various parameters such as temperature and chemical 
potential. 



2 Data set and model description 

As emphasized in the introduction, in the present analysis we use the most recent available data, 
concentrating on fully integrated particle yields and discarding data that have been obtained 
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in limited kinematic windows. The only exceptions to this rule are the A/ A and p/K~ ratios 
in Si-Au collisions |12, [lj| which were not available in full phase space. It has been decided to 
keep them as they are the only available recent measurements involving antibaryons. 

We have derived integrated multiplicities of tt + , A and proton in Au-Au collisions at AGS by 
extrapolating published rapidity distributions [1IBL O, [TBI with constrained mid-rapidity value 



(z/nn=1-6). For proton and A we have fitted the data to Gaussian distributions, whilst for n + 
we have used a symmetric flat distribution at midrapidity with Gaussian-shaped wings on each 
side; the point at which the Gaussian wing and the plateau connect is a free parameter of the fit. 
The fits yielded very good x 2 's/dof: 0.27, 1.24 and 1.00 for tt + , proton and A respectively. The 
integrated multiplicities have been taken as the area under the fitted distribution between the 
minimal y m i n and maximal y max values of rapidities for the reactions NN — > 7rNN, NN — > AK for 
pions and A's respectively; the difference between these areas and the total area has been taken 
as an additional systematic error. The area between y min and y max amounts to practically 100% 
of the total area for pions and about 95% for A's. Ref. [IT] quotes an additional experimental 
systematic error of 10% on A multiplicity that we have added in quadrature. Hence we obtain: 



(tt + ) = 133.7 ±9.93 
(A) = 20.34 ±1.36 ±1.23 ±2.03 (1) 

where the first error is the fit error, the second is the systematic error due to the variation 
of integration region and the third is the experimental systematic error. As to protons, the 
extracted rapidity interval corresponding to the reaction N N — > N N is only 79% of the 
total Gaussian area. The difference between the two areas is too large to be considered as an 
additional error; thus, in order to reduce the uncertainty, we have decided to take the ratio 
p/7r + extracted in the above rapidity interval rather than the proton multiplicity itself. This 
yields: 

(p/tt+) = 1.234 ±0.126 (2) 

where the error includes both the fit error and an error stemming from a 10% systematic 
uncertainty quoted in ref. WE . 



We have not included data on deuteron production because of the possible inclusion of 
fragments in the measured yields. This is particularly dangerous at low (SIS) energies where 
inclusion or not of deuterons modifies thermodynamic quantities like e/n ||19|| . 

The data analysis has been performed within an ideal hadron gas grand-canonical framework 
for Pb-Pb at SPS and Au-Au at AGS whereas for Au-Au at SIS and Si-Au at AGS we have 
required the exact conservation of strangeness instead of using a strangeness chemical potential 
(see the discussion later in the text); in both cases we have used a supplementary strange quark 
fugacity 7 S . In the grand-canonical approach, the overall average multiplicities of hadrons and 
hadronic resonances are determined by an integral over a statistical distribution: 

(m) = (2J, + 1)7^5 / d 3 p — -— i — — — (3) 

(2ny J 7 S s * exp [(£* - ^ • q.i)/T\ ± 1 

where is a three-dimensional vector with electric charge, baryon number and strangeness of 
hadron i as components; fi the vector of relevant chemical potentials; J{ the spin of hadron i 
and Si the number of valence strange quarks in it; the ± sign in the denominator is relevant 
for fermions, the — for bosons. This formula holds in case of many different statistical-thermal 
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systems (i.e. clusters or fireballs) having common temperature and 7 S but different arbitrary 
momenta, provided that the probability of realizing a given distribution of quantum numbers 
among them follows a statistical rule [0, |2"0fl . In this case V must be understood as the sum 
of all cluster volumes measured in their own rest frame. Furthermore, since both volume and 
participant nucleons may fluctuate on an event by event basis, V and n (and maybe T) in 
Eq. (|) should be considered as average quantities 0. 

The overall abundance of a hadron of type % to be compared with experimental data is 
determined by the sum of Eq. (|3|) and the contribution from decays of heavier hadrons and 
resonances: 



primary 

n,- = nf 



+ ]T Br(j -> i) nj (4) 
j 

where the branching ratios Br(j — > i) have been taken from the 1998 issue of the Particle Data 



Table 21 



It must be stressed that the unstable hadrons contributing to the sum in Eq. (|4]) may differ 
according to the particular experimental definition. This is a major point in the analysis pro- 
cedure because quoted experimental multiplicities may or may not include contributions from 
weak decays of hyperons and K^. We have included all weak decay products in our computed 
multiplicities except in Pb-Pb collisions on the basis of relevant statements in ref. [02] and 



about antiproton production in refs. ]15[ It must be noted that switching this assumption 
in Au-Au at SIS and AGS does not affect significantly the resulting fit parameters whereas it 
does in Si-Au. 

The overall multiplicities of hadrons depend on several unknown parameters (see Eq. (|3|)) 
which are determined by a fit to the data. The free parameters in the fit are T, V, 7 S and \ib 
(the baryon chemical potential) whereas fj,s and fiQ, i.e. the strangeness and electric chemical 
potentials, are determined by using the constraint of overall vanishing strangeness and forcing 
the ratio between net electric charge and net baryon number Q/B to be equal to the ratio 
between participant protons and nucleons. The latter is assumed to be Z/A of the colliding 
nucleus in Au-Au and Pb-Pb while it has been calculated to be 0.43 for central Si-Au collisions 
by means of a geometrical model. 

As we have mentioned before, for SIS Au-Au and AGS Si-Au data we have required the 
exact conservation of strangeness instead of using a strangeness chemical potential. This gives 
rise to slightly more complex calculations |24| which are necessary owing to either very small 
strange particle production (Au-Au) or a relatively small system size (Si-Au). The difference 
between the strangeness-canonical and pure grand-canonical calculations of multiplicities of K 
and A for the final set of thermal parameters (see Table 1) turns out to be around 2-3% for K 
and A in Si-Au but it is as large as a factor 15 in Au-Au at 1.7 A GeV/c. 

Owing to few available data points in SIS Au-Au collisions, we have not fitted the volume V 
nor the 7 S therein. The volume has been assumed to be 47rr 3 /3 where r = 7 fm (approximately 
the radius of a Au nucleus) while 7 S has been set to 1, the expected value for a completely 
equilibrated hadron gas. Since we have performed a strangeness-canonical calculation here, 
the yield ratios involving strange particles are not independent of the chosen volume value as 
in the grand-canonical framework. Thus, in this particular case, V is meant to be the volume 
within which strangeness is conserved (i.e. vanishing) and not the global volume defining overall 
particle multiplicities as in Eq. (Q). Also, in order to test the dependence of this assumption 
on our results, we have repeated the fit by varying V by a factor 2 and 0.5 in turn. 

The yields of resonances have been calculated by integrating Eq. (§) times a relativis- 
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tic Breit-Wigner distribution over an interval [m — 8mi,m + 5m u ], where 8mi = min[m — 
^threshold, 2 r] and 5m u = 2T. The minimum mass m threshold is required to open all decay 
modesQ The relativistic Breit-Wigner distribution has been renormalised within the integra- 
tion interval. The non-vanishing width of resonances plays a major role especially at low 
energies (e.g. SIS); for instance, the A (1232) resonance creates pions more effectively than in 
the case of a vanishing width. 

We have not used proper volume corrections in a Van der Waals type fashion which have 
been considered previously p5| . 

A major problem in Eq. (Q) is where to stop the summation over hadronic states. Indeed, 
as mass increases, our knowledge of the hadronic spectrum becomes less accurate; starting 
from rs 1.7 GeV many states are possibly missing, masses and widths are not well determined 
and so are the branching ratios. For this reason, it is unavoidable that a cut-off on hadronic 
states be introduced in Eq. (|). If the calculations are sensitive to the value of this cut-off, 
then the reliability of results is questionable. We have performed all our calculations with two 
cut-offs, one at 1.8 GeV (in the analysis algorithm A) and the other one at 2.4 GeV (in the 
analysis algorithm B). The contribution of missing heavy resonances is expected to be very 
important for temperatures > 200 MeV making thermal models inherently unreliable above 
this temperature. 



3 Results 

As mentioned in the introduction, we have performed two analyses (A and B) by using com- 
pletely independent algorithms. In the analysis A all light- flavoured resonances up to 1.8 GeV 
have been included. The production of neutral hadrons with a fraction / of ss content has 
been suppressed by a factor (1 — /) + fj^. In the analysis B the mass cut-off has been pushed 
to 2.4 GeV and neutral hadrons with a fraction / of ss content have been suppressed by a 
factor jg* . Both algorithms use masses, widths and branching ratios of hadrons taken from the 
1998 issue of Particle Data Table [plfl . However, it must be noted that differences between the 
two analyses exist in dealing with poorly known heavy resonance parameters, such as assumed 
central values of mass and width, where the Particle Data Table itself gives only a rough esti- 
mate. Moreover, the two analyses differ by the treatment of mass windows within which the 
relativistic Breit-Wigner distribution is integrated. 

The results of the \ 2 fits are shown in Tables 1 and 2 for both analyses A and B. The 
agreement is indeed very good and confirms the reliability of the results obtained. The x 2 
minimisation in Au-Au collisions at AGS in analysis B did not converge to a reliable minimum; 
however, the \ 2 computed in analysis B by fixing the values of the parameters to the ones 
obtained in analysis A is approximately as large as in analysis A itself, thus confirming the 
good agreement between the two calculations. 

We have investigated in detail the lack of convergence of analysis B in Au-Au collisions. The 
main reason of the fragility of the fit is the absence of measured antibaryon yields, which are 
very effective in fixing the baryon-chemical potential, in the main set of full phase space data. 
That shortage brings about a shallowness of \ 2 minima in four dimensions, and, consequently, 
a difficult convergence in both analyses. Notwithstanding, in analysis A the absolute minimum 
turned out to be deep enough, whereas in analysis B the convergence to a sufficiently nearby 
point was spoiled and the minimum drifted to T « 150 MeV with a nearly flat descent from 



fact, in analysis A (see below) the integration interval has been taken symmetric [m — 8mi,m + 5rni 
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the minimum found in A. This indicated a possible model dependence of the fit outcome. In 
order to check our result in analysis A and make it robust we have repeated the fits in Au-Au 
collisions at AGS by using an additional measurement of p/p|23| ratio in the limited phase space 
region 1.0 < y < 2.2 around midrapidity. The use of a ratio of particles measured by the same 
experiment under the same conditions reduces the involved systematic errors due to slightly 
different centrality definitions (with respect to the other data set) and other possible sources. 
However, the actual ratio in full phase space might be different owing to different shapes of 
p and p rapidity distribution and this effect has been taken into account by conservatively 
assigning a 20% additional systematic error. The fit results are shown in Table 3; the two 
analyses are in very good agreement and, on top of that, the results for analysis A are in 
excellent agreement with those in Table 1 obtained without using p/p ratio, thus confirming 
the good quality of the calculation. 

For each analysis an estimate of systematic errors on fit parameters have been obtained by 
repeating the fit 

• assuming vanishing widths for all resonances 

• varying the mass cut-off to 1.7 GeV in analysis A and to 1.8 GeV in analysis B 

• for Au-Au at 1.7 A GeV/c, the volume V has been varied to V/2 and to 2V (see discussion 



The differences between the new fit parameters and the main parameters have been conser- 
vatively taken as uncorrelated systematic errors to be added in quadrature for each variation 
(see Table 1). The effect of errors on masses, widths and branching ratios of inserted hadrons 
has been studied in analysis A according to the procedure described in ref. |7j] and found to be 
negligible. 

Finally, the results of the two analyses have been averaged according to a method suggested 
in ref. |3"0|] , well suited for strongly correlated measurements. Firstly, a simple no-correlation 
weighted average has been calculated as the central value of each parameter. Secondly, the 
error on it has been estimated by conservatively assuming that the results A and B are fully 
correlated, i.e. with a covariance matrix: 



The correlation between analyses A and B clearly arises from the use of the same set of hadronic 
data and theoretical model. 

In Table 1 we also list the values of the Wroblewski factor A s [|3l| measuring the number of 
newly created primary valence ss pairs in comparison to the newly created non-strange primary 
valence quark pairs 



in Sect. 2) 




(5) 



yielding an error: 




(6) 



2(ss) 



(7) 
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along with fit and systematic errors. The ss and light quark pairs are computed on the basis 
of primary multiplicities of all hadron species, i.e. before particle decays take place. The 
behaviour of A s as a function of collision type and centre-of-mass energy is shown in Fig. |1] 
including elementary and S-S, S-Ag collisions. Values for S-S, S-Ag and e + e~, pp, pp collisions 
have been taken from ref. [^]. 



4 Discussion and conclusions 

From the results obtained, it emerges that a statistical-thermal description of multiplicities in 
a wide range of heavy ion collisions is indeed possible to a satisfactory degree of accuracy, for 
beam momenta ranging from 1.7 A GeV/c to 158 A GeV/c per nucleon. Furthermore, the fitted 
parameters show a remarkably smooth and consistent dependence as a function of centre-of- 
mass energy. The fit quality is generally good with the exception of Au-Au collisions at 1.7 A 
GeV/c where the large x 2 is due to an underestimation of one ratio tj/tt (see Table 2). 

The temperature varies considerably between the lowest and the highest beam energy, 
namely, between 50 MeV at SIS and 160 MeV at SPS. Similarly, the baryon chemical po- 
tential changes appreciably, decreasing from about 820 MeV at SIS to about 240 MeV at SPS. 
However, since the changes in temperature and chemical potential are opposite, the resulting 
energy per particle shows little variation and remains practically constant at about 1 GeV per 
particle; this is shown in Fig. ^|. 

The supplementary 7 S factor, measuring the deviation from a completely equilibrated hadron 
gas, is around 0.7 - 0.8 at all energies where it has been considered a free fit parameter. At 
the presently found level of accuracy, a fully equilibrated hadron gas (i.e. 7 S =1) cannot be 
ruled out in all examined collisions except in Pb-Pb, where 7 S deviates from 1 by more than 
4a. This result does not agree with a recent similar analysis of Pb-Pb data Jl3|| imposing a full 



strangeness equilibrium. The main reason of this discrepancy is to be found in the different 



data set used; whilst in ref. [13] measurements in different limited rapidity intervals have been 
collected, we have used only particle yields extrapolated to full phase space. The temperature 
values that we have found essentially agree with previous analyses in Au-Au collisions |26fl and 



Si-Au collisions ||T0| , [13| , p2| and estimates in Au-Au collision at 11.7 A GeV/c ||33|| . 

The T value in Pb-Pb is significantly affected by the multiplicity value of the heaviest 
particles measured, namely <fi and 5, as they are almost entirely directly produced and provide 
a major lever arm on the slope of production vs. mass function. A recent 40% lowering of 
the S~ yield measured by NA49 ]29| with respect to a previous measurement [34] results in a 



decrease of estimated temperature value from about 180 MeV to the actual 160 MeV. However, 
the removal of these two particles from the data set yields fitted parameter values which are 
in fair agreement with the main fit, as shown in Table 3. In particular, it is worth remarking 
that this exclusion does not bring significant change to 7 S whose outcome is very sensitive to 
particles with multiple strange quark content and this confirms again the robustness of the 
main fit. 

In order to further investigate strangeness production in Pb-Pb we have performed a consis- 
tency test between our fitted parameters, based on NA49 measurements, and the multiplicities 
of multi-strange hadrons measured by the experiment WA97 in central Pb-Pb collisions in a ra- 
pidity window ±0.5 around mid-rapidity []3"5| . By fixing T, 7 S and /ie to the averaged values in 
Table 1 and adjusting the volume (i.e. an overall normalisation), we obtain a x 2 /dof = 28.9/6. 
Calculated A, A multiplicities (see Table 4) do not include a residual feeding from S decays in 
the experiment, estimated to be < 5%, < 10% respectively ^6[. The high value of the \ 2 indi- 
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cates that the statistical-thermal analysis is not able to reproduce data in a limited phase space 
region and in full phase space at the same time without resorting to a more detailed dynamical 
model. In particular, the parameters determined by the fit to NA49 data underestimate the 
yields of 5 and Q baryons. 

The parameter 7 S as a function of centre-of-mass energy in heavy ion collision (including S-S 
and S-Ag 0) is shown in Fig. |^. Again, the values for S-S, S-Ag and e + e~, pp, pp collisions 
have been taken from ref. J7|. As can be seen from the Fig. 7 S is fairly constant, however, 
given the large error bars, it is quite difficult to exclude different behaviours. Also the behaviour 
of A s factor (see Fig. [I]) as a function of energy (provided that there is little dependence on 
system size at fixed i/i, as the approximate equality of A s in S-S and S-Ag confirms) is still 
unclear due to large experimental uncertainties. The line shape is either compatible with a 
monotonically increasing curve, saturating at A s ~ 0.45, or with a curve having a maximum 
around Si-Au collisions, then decreasing and settling at an asymptotic ~ 0.45 value or maybe 
decreasing further to the characteristic value of elementary collisions. 

Forthcoming lower energy Pb-Pb and high energy Au-Au data at RHIC should allow to 
clarify the behaviour of strangeness production in heavy ion collision. In order to easily com- 
pare our results with new measurements from RHIC experiments we also show in figure 4 the 
values of various particle-antiparticle ratios as a function of p/p ratio for different values of 
the temperature (T=160, 165 and 170 MeV) and a fixed value of the charge to baryon ratio of 
0.401. RHIC results, however, will only be available for very limited kinematical region, while 
this kind of thermal model approach is largerly tied to full phase space ratios (see Introduction). 
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Table 1: Summary of fit results. Free fit parameters are quoted along with resulting minimum 
X 2 's and A s parameters. 





Analysis A 


Analysis B 


Average 


Au-Au 1.7 A GeV 


T (MeV) 


49.6±1.0±2.2 


49.7±1.1±2.3 


49.6±2.5 


Hb (MeV) 


810±15±12 


818±15±12 


813±23 


7s 
\/(fm 3 ) 


1 (fixed) 
1437 (fixed) 


1 (fixed) 
1437 (fixed) 


1 (fixed) 
1437 (fixed) 


X Vdof 


14.9/2 


15.1/2 




As 


0.0050±0.0034 


0.0058±0.0036 


0.0054±0.0035 


Au-Au 11.6,4 


GeV 


T (MeV) 


121.2±3.9±3.0 




121.2±4.9 


\ib (MeV) 


559±15±9 




559.4±16 




0.697±0.080±0.043 




0.697±0.091 


\/T 3 exp(-0.7GeV/T) 


2.01±0.23±0.14 




2.01±0.27 


X 7dof 
As 


2.25/2 
0.43±0.10 




0.43±0.10 


Si-Au 14.64 GeV 


T (MeV) 


133.1±3.9±1.4 


138.0±4.5±0.3 


135.4±4.3 


Hb (MeV) 


592±34±13 


573±30±0.8 


581±32 




0.843±0.095±0.067 


0.847±0.061±0.069 


0.845±0.101 


\/T 3 exp(-0.7GeV/T) 


0.526±0.090±0.081 


0.545±0.079±0.12 


0.534±0.130 


X 2 /dof 


14.3/4 


11.6/4 




As 


0.74±0.2 


0.72±0.12 


0.72±0.14 


Pb-Pb 158 A GeV 


T (MeV) 


159.5±2.5 ±1.5 


156.0±2.4±2.6 


158.1±3.2 


/xb (MeV) 


238 ±13 ±3 


239±12±5 


238±13 


7s 


0.760±0.035±0.028 


0.862±0.036±0.061 


0.789±0.052 


VT 3 exp (-0.7GeV/T) 


20.9 ±1.5 ±2.0 


19.7±1.0±2.9 


21.7±2.6 


X 2 /dof 
As 


14.4/6 
0.444±0.026 


22.6/6 
0.450±0.024 


0.447±0.025 
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Table 2: Comparison between fitted and measured particle multiplicities and ratios. In Au-Au 
collisions at AGS we also quote our prediction (including weak decay products) along with a 
measurement in a limited kinematic range (1.0 < y < 2.2), whose error is only statistical, which 
sets a lower limit for p multiplicity. 





Reference 


Measurement 


Analysis A 


Analysis B 


Au-Au 1.7 A GeV 


7T+/p 






0.052±0.013 


0.05306 


0.05306 


K+/V+ 




2§ 


0.003±0.00075 


0.003040 


0.003030 


7r~ /ir + 




2§ 


2.05±0.51 


2.0371 


2.007 


ri/ir° 




2§ 


0.018±0.007 


0.00109 


0.000851 


Au-Au 11. 6 A GeV 


Participants 






363±10 


363.0 




K+ 






23.7±2.9 


20.23 




K~ 






3.76±0.47 


4.038 




7T+ 




nj 


133.7±9.9 


133.3 




A 




ij 


20.34±2.74 


21.54 




p/vr+ 


1 


1 


1.234±0.126 


1.295 




P 




1 


>0.0185±0.0018 


0.0363 




Si-Au 14.64 GeV 


Participants 






115±10 


99.41 


94.33 


7T+ 




2§ 


33±3 


34.79 


36.92 


7r~ /ir + 




2§ 


1.09±0.13 


1.296 


1.196 


K+/7T+ 




if 


0.18±0.02 


0.1564 


0.1590 


K _ /ir _ 




3 


0.034±0.004 


0.02715 


0.02767 


p/K- 






0.018±0.0034 


0.01672 


0.01706 


A/A 




1 


0.003±0.0015 


0.00217 


0.00301 








0.09±0.04 


0.1725 


0.1345 




Pb-Pb 158A GeV 


(7r+ + 7r-)/2. 






600±30 


581.9 


568.0 


K+ 




H 


95 ±10 


96.42 


99.05 


K" 




3 


50 ±5 


56.53 


60.96 






21 


60 ±12 


75.39 


79.34 


P 




2l 


140±12 


144.9 


144.9 


P 




11 


10 ±1.7 


8.242 


7.707 


„ 






7.6±1.1 


7.185 


5.852 






2§ 


4.42±0.31 


3.895 


4.110 






n 


0.74±0.04 


0.766 


0.765 


A/A 




2| 


0.2±0.04 


0.1033 


0.098 
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Table 3: Top: fit results for Au-Au collisions with the measured ratio p/p; its systematic error 
has been conservatively estimated to be 20%. Bottom: fit results for Pb-Pb collisions with the 
exclusion of 0, S's and both. 



Au-Au 11.64 GeV 





Analysis A 


Analysis B 


Measurement 


T (MeV) 


121.2±4.6±1.7 


130.6±5.5±3.9 




He (MeV) 


558±12±9 


594±26±30 




7s 


0.701±0.068±0.072 


0.883±0.124±0.207 




VT 3 exp (-0.7GeV/T) 


2.02±0.19±0.36 


1.65±0.22±0.66 




X 2 /dof 


2.70/3 


1.06/3 




As 


0.44±0.11 


0.72±0.30 




Participants 


363.0 


364.1 


363±10 


K+ 


20.23 


21.37 


23.7±2.9 


K" 


4.038 


3.950 


3.76±0.47 


7T + 


133.3 


130.5 


133.7±9.9 


A 


21.54 


21.40 


20.34±2.74 


p/vr + 


1.295 


1.237 


1.234±0.126 


p/p 


2.125 10~ 4 


2.516 10" 4 


(2.50±0.25±0.50)10* 4 


Pb-Pb 1584 GeV - Analysis A 




without 


without S's 


without and S's 


T (MeV) 


159.9±2.6 


168.7±7.1 


165.6±6.5 


fJ>B (MeV) 


237 ±13 


222±15 


218±9 


Ts 


0.753±0.036 


0.710±0.051 


0.664±0.064 


VT 3 exp (-0.7GeV/T) 


22.8 ±1.2 


23.8±1.3 


24.4±1.3 


X 2 /dof 


14.0/5 


5.5/4 


4.5/3 





7.092 


6.829 


6.053 




3.849 


3.616 


3.154 




0.7637 


0.9749 


0.8028 


Pb-Pb 1584 GeV - Analysis B 




without 


without S's 


without and S's 


T (MeV) 


158.2±2.7 


167.1±6.3 


158.5±5.2 


He (MeV) 


232±12 


227±13 


208±14 




0.806±0.040 


0.862±0.043 


0.658±0.067 


VT 3 exp (-0.7GeV/T) 


20.9 ±1.1 


20.2±1.1 


22.9±1.3 


X 2 /dof 


15.9/5 


14.6/4 


6.3/3 




4.784 


6.293 


2.328 




3.867 


4.397 


2.661 




0.7594 


1.151 


0.5916 
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Table 4: Comparison between predicted particle multiplicities in central Pb-Pb collisions by 
using fitted T, /ib and 7 S quoted in the rightmost column of Table 1 and those measured by 
WA97. The normalisation volume has been adjusted to minimise a x 2 which turned out to be 
28.9. 





Calculated 


Measured 


h- 

A 
A 

n + n 


207.9 
23.66 
15.58 
1.543 
1.251 
0.2354 
0.1662 


178±22 
21.9±2.4 
13.7±0.9 
1.8±0.2 
1.5±0.1 
0.37±0.06 
0.41±0.08 
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Figure 1: Ratio of strange quark pairs to created non-strange quark pairs (Wroblewski factor) 
as a function of nucleon-nucleon centre-of-mass energy. Values for S-S, S-Ag and e + e~, pp, 
pp collisions have been taken from ref . . 
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Figure 2: Fitted temperatures and baryon-chemical potentials plotted along with curves of 
constant energy per hadron. 
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Figure 3: Strangeness suppression factor 7 S as a function of nucleon-nucleon centre-of-mass 
energy. Values for S-S, S-Ag and e + e~, pp, pp collisions have been taken from ref. M 
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Figure 4: Particle-Antiparticle and negative to charged hadrons ratios as a function of the p/p 
ratio for different temperatures and a fixed ratio of charge over baryon number Q/B 0.401 and 
7 S =1. The p/p variation is governed by the variation of the baryon chemical potential. The 
dependence of these ratios on 7 S , within the expected range, is found to be negligible. 
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